Heterogeneous aging in spin glasses 
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We introduce a set of theoretical ideas that form the basis for an analytical framework capable of 
describing nonequilibrium dynamics in glassy systems. We test the resulting scenario by comparing 
its predictions with numerical simulations of short-range spin glasses. Local fluctuations and re- 
sponses are shown to be connected by a generalized local out-of-equilibrium fluctuation-dissipation 
relation. Scaling relationships are uncovered for the slow evolution of heterogeneities at all time 
scales. 
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PACS: 75.10.Nr, 75.10.Jm, 75.50.Lk, 75.10.Hk 

Very slow equilibration and sluggish dynamics are 
characteristics shared by disordered spin systems and by 
other glassy systems such as structural and polymeric 
glasses. The origin of this dynamic arrest near and be- 
low the glass transition is currently poorly understood. 
Studies of the time evolution of many quantities, such 
as the remanent magnetization, the dielectric constant, 
or the incoherent correlation function, have shown that 
below the glass transition the system falls out of equilib- 
rium This is evidenced by the presence of aging, 
i.e. the dependence of physical properties on the time 
since the quench into the glassy state, and also by the 
breakdown of the equilibrium relations dictated by the 
fluctuation dissipation theorem (fdt) 

Most analytical progress in understanding non- 
equilibrium glassy dynamics has been achieved in mean- 
field fully connected spin models g], while numerical 
simulations have addressed both structural glasses ||] 
and short-range spin glass models Until recently, 
however, experimental, numerical, and analytical studies 
have mainly focused on global quantities, such as global 
correlations and responses, which do not directly probe 
local relaxation mechanisms. Local regions that behave 
differently from the bulk, or dynamic heterogeneities, 
could be crucial to understand the full temporal evolu- 
tion, and have received considerable experimental 
and numerical |^p[] attention lately. However, no clear 
theoretical picture has yet emerged to describe the local 
nonequilibrium dynamics of the glassy phase. 

Here we introduce such a theoretical framework, and 
test its predictions via numerical simulations of a short- 
range spin glass model. We show that local correlations 
and responses are linked, and we find scaling properties 
for the heterogeneities that connect the evolution of the 
system at different times. This universality may provide 
a general basis for a realistic physical understanding of 
glassy dynamics in a wide range of systems. 

The framework that we propose is motivated by an 
analogy ]TT|] between aging dynamics and the well-known 
statics of Heisenberg magnets. For concreteness, we 



test its predictions against Monte Carlo simulations on 
the prototypical spin glass model, the three dimensional 
Edwards-Anderson (3DEA) model, H = Yl{ij) JijSiSj, 



where 



±1 and the nearest-neighbor couplings are 



Jij = ±1 with equal probability. We argue that two 
dynamical local quantities, the coarse-grained local cor- 
relation C r (t,t w ) = Y^2 ieV Si(t)si(t w ) and integrated 

response t w ) = -^j J^kli V ^ieV r "'^w — are 
essential to understand the mechanisms controlling the 
dynamics of glassy systems. The spins are represented 
by Si in the absence of an applied field and by Si\ h (k) in 

the presence of one. ~§i(t) = ^EtyV'f'') i s tne re ~ 
suit of coarse-graining the spin over a small time-window 
[typically, r = 1000 Monte Carlo steps (MCs)]. V r is a 
cubic box with volume V centered at the point r. By 
taking V to be the volume of the whole system, the bulk 
or global correlation C(t,t w ) and response x(t, t w ) are 
recovered. Two generic times after preparation are rep- 
resented by t w and t, with t w < t. When the system is 
not in equilibrium, time dependences do not reduce to a 
dependence on the time difference t — t w . We measure a 
staggered local integrated linear response by applying a 

(k) 

bimodal random field on each site h\ ' — ±h during the 
time interval [t w , t] . Linear response holds for the values 
of h that we use. The index k = 1, . . . , Nt labels different 
realizations of the perturbing field. We use random ini- 
tial conditions. The thermal histories, i.e. the sequences 
of spins and random numbers used in the MC test, are 
the same with and without a perturbing field. 

In a disordered spin model, the coarse-grained local 
magnetization typically vanishes, but the local corre- 
lation is non-trivial. Averaged over disorder and the 
thermal history, this correlator defines the Edwards- 
Anderson parameter q EA when t w — > oo, and t — t w —> oo 
subsequently. Can we detect the growth of local or- 
der by analyzing the evolution of the local correlator, 
as one easily can for a system undergoing ferromagnetic 
domain growth? In Fig. |l| we show the local correlation 
for fixed t w and t on a 2D cut of the 3DEA model. Re- 
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gions with large values of C r are intertwined with regions 
with a small value of C r as shown by the contour levels. 
This behavior persists for all t w and t that we can reach 
with the simulation and a more sophisticated analysis is 
necessary to identify a growing order in this system. 




FIG. 1. The local correlations on a 2D cut of the cubic cell. 
The linear size is L = 64. T = 0.72 T c , V = 3 3 , t w = 4.1 x 10 5 
MCs and t = 2.8 x 10 6 MCs. 

It is clear from Fig. [l] that different sites have distinct 
dynamics. Analysis of the local correlation for fixed t w 
as a function of t shows that in general the relaxation is 
non-exponential - this is often ascribed to the presence of 
heterogeneous dynamics. How can one characterize the 
heterogeneous dynamics and determine its origin? We 
argue that relevant quantities are the probability distri- 
bution (pdf) of the local correlation, p(C r (t,t w )), the 
PDF of the local integrated response, p(x r (t,t w )), and 
the joint PDF p(C r (t,t w ),Xr(t,t w )) and we start by dis- 
cussing the latter. 

The fdt relates the correlation of spontaneous fluctua- 
tions to the integrated linear response of a chosen observ- 
able (e.g. C r (t,t w ) and Xr(t,t w ) averaged over thermal 
histories) at equilibrium. Glassy systems modify the fdt 
in a particular way first obtained analytically for mean- 
field models , later verified numerically in a number of 
realistic models |l^,Q, and more recently tested exper- 
imentally [p| . A parametric plot of the bulk integrated 
response, x(t, t w ), against the bulk correlation, C(t,t w ), 
for fixed and long waiting-time t w and using t as a param- 
eter, approaches a non-trivial limit, x(C), represented by 
the crosses in Fig. ||(b). The curve has a straight section, 
for which t — t w <^Lt w and the correlation decays from 1 
to q EA with a slope — 1/T as found with the equilibrium 
fdt. Beyond this point, as t increases towards infin- 
ity, the curve separates from the fdt line. Now, consider 
each lattice site for fixed times t w and t. If we plot points 
for the pairs (C r (t,t w ),Xr(t,t w )), where will they lie? 

When t w — > oo and t — t w -C t w> all local correlators 
satisfy the fdt strictly once averaged over thermal histo- 
ries, since the magnitude of local deviations from the fdt 
has an upper bound |fl6|| . We have checked that C r (t, t w ) 



and Xr(t, tw) obey the fdt for an individual thermal his- 
tory apart from small fluctuations [see Fig. §(b)]. 

For the regime of widely separated times we propose an 
analysis similar in spirit to the one that applies to the low 
energy excitations of the Heisenberg model. There, the 
free energy for the coarse grained magnetization m(r) is 
F = J d d r[{Vprfi(r)) 2 + T^(|m(r)|)]. A spontaneous sym- 
metry breaking signals the transition into the ordered 
phase (m) = mo ^ 0, in which the order parameter 
has both a uniform length (the radius of the bottom of 
the effective potential V(|m|)), and a uniform direction. 
F is invariant under uniform rotations rh(r) —* lZm(r). 
The lowest energy excitations (spin waves) are obtained 
from the ground state by leaving the length of the vec- 
tor invariant and applying a slowly varying rotation to 
it: rh(r) = TZ(r)mQ. These are massless transverse fluc- 
tuations (Goldstone modes). In contrast, longitudinal 
fluctuations, which change the magnitude of the magne- 
tization vector, are massive and energetically costly. 

Let us now apply the same kind of analysis to the dy- 
namics of the spin glass. Here, the relevant fluctuating 
quantities are the coarse grained local correlations C r and 
their associated local integrated responses Xr ■ In Rcf. Jll]] 
we derived an effective action for these functions that 
becomes invariant under a global time-reparametrization 
t — *■ h{t) in the aging regime. This symmetry leaves the 
bulk relation, x(C), invariant. A uniform reparametriza- 
tion is analogous to a global rotation in the Heisenberg 
magnet, and the curve x(C) is analogous to the surface 
where V(|m|) is minimized. Hence, we expect that for 
fixed long times t w and t in the aging regime, the lo- 
cal fluctuations in C' r and Xr should be given by smooth 
spatial variations in the time reparametrization, h r (t), 
i.e. C r (t,t w ) = C SP (h r (t),h r (t w )) « C(h r (t),h r (t w )) 
where C SP is the global correlation at the saddle-point 
level that in the numerical studies we approximate by 
the actual global correlation C, and similarly for Xr- 
These transverse fluctuations are soft Goldstone modes. 
Longitudinal fluctuations, which move away from the 
x(C) curve, are massive and penalized. This implies the 
first testable prediction of our theoretical framework: the 
pairs (C r ,Xr) should follow the curve x(C) f° r the bulk 
integrated response against the bulk correlation. 

In Fig. ||, we test this prediction by plotting the dis- 
tribution of pairs (C r ,Xr)- We find, as expected, that 
for long times the dispersion in the longitudinal direc- 
tion (i.e. away from the bulk x(C) curve) is much weaker 
than in the transverse direction (i.e. along the bulk x(C) 
curve). In the coarse grained aging limit we expect the 
former to disappear while the latter should remain. (This 
limit corresponds to the way actual measurements are 
performed: the thermodynamic limit is taken first to 
eliminate finite size effects and undesired equilibration; 
then the large t w limit is taken to reach the asymptotic 
regime; finally, the limit V — > oo serves to eliminate fluc- 
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tuations through the coarse graining process; in the figure 
we used a large volume V — 13 3 to approach the latter 
limit though we found a similar qualitative behavior for 
smaller V.) Figure ||(a) displays the joint PDF p(C r ,Xr) 
for a pair of times (t w , t) that are far away from each 
other. Figure 0(b) shows the projection of a set of con- 
tour levels for t w fixed and six values of t. Even though 
the data for each contour corresponds to a single pair 
of times (t w , t) , the fluctuations span a range of values 
that, for the bulk quantities, would require a whole fam- 
ily of pairs (t w , t) . This reveals that the aging process is 
non-uniform across a finite-range model. 
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FIG. 2. L = 64, T = 0.72 T c , V = 13 3 , t w = 4 X 10 4 . (a) 
Surface plot of the joint PDF for t/t w = 16. Points within the 
shown contour account for 2/3 of the total probability, (b) 
Time evolution of p(C r (t, t w ), Xr(t, in,)), for t/t w = 1.00005, 
1.001, 1.06, 2, 8, 32 (from right to left). Each contour contains 
2/3 of the total probability for the given pair of times. The 
line marked fdt is the equilibrium relation between \r and 
C r - The crosses are the bulk x(C) f° r the same pairs of times. 

We now turn our attention to a more detailed exami- 
nation of the time dependences in the local dynamics. A 
good t/t w scaling that breaks down only for very large 
values of the subsequent time t — t w has been obtained 
for the bulk thermoremanent magnetization experimen- 
tally j|] and for the bulk correlation numerically || once 
the stationary (t — t w <C t w ) part of the relaxation is 



subtracted, as suggested by the solution to mean-field 
models M. For systems that display this particular de- 
pendence on t/t w for the bulk correlator, a second pre- 
diction can be extracted from our theoretical framework: 
the distribution p(C r (t,t w )) should only depend on the 
ratio t/t w . Even further, if the bulk correlator has a 
simple power law form C SP (t,t w ) ~ q E A(t/t w )~ p , an ap- 
proximate treatment of fluctuations leads to a rescaling 
and collapse of p(C r (t,t w )) even for pairs of times with 
different ratios t/t w . 

Since we are dealing with ratios of times, it is con- 
venient to define h r (t) = e Vr ^ t \ so that C r (t,t w ) = 
C SP (h r (t)/h r (t w )) = ^(e^O-^ 1 -)). Therefore the 
statistics of local correlations are determined from the 
statistical distribution of distances between two "sur- 
faces", <p r (t) — ip r (t w ). In this form, a dynamic theory 
of short-range spin glasses is not different from a the- 
ory of fluctuating geometries or elasticity. We propose 
a simple reparametrization invariant effective action for 
(fi r (t) = hit + Sip r (t), expanding around 5ip r (t) — 0, with 
no zeroth or first order term in 6tp r (t). We assure that the 
effective action is reparametrization invariant by taking 
one time derivative for each time variable. Thus fllTtl 
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d d r / dt / dt' X7tp r (t) V(p r (t') e -p\vAt)-Vr{t')\^ 



where the last factor penalizes fast time variations of f r 
and the V ensure that spatial variations are smooth. 
Expanding to lowest order in 5(p r (t) yields (p r (t) — 
<Pr(tw) = ln(t/t w ) + 6<p r (t) - 5(p r (t w ) ~ \n(t/t w ) + (a + 
b\n(t/t w )) a X r (t,t w ), where a and b are determined by 
the magnitude of the fluctuations, and X r (t, t w ) is a ran- 
dom variable drawn from a time-independent PDF that 
governs the fluctuations of the surfaces. In our approx- 
imation, which describes uncorrelated drift between two 
surfaces (i.e. a random walk), a = 1/2. 
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FIG. 3. PDF of the local correlations for the ratios 
t/t w = 2, 8, and 32, averaged over 36 realizations of disor- 
der. L = 32, T = 0.72 T c and V = 3 3 . The waiting-times are 
t w = 10 4 (2™™" 1 ) MCs, with n w given in the key. We include 
the PDF of the local integrated responses, Xr(t, t w ), for n w — 3 
and t/tw = 2. (p[Xr) is divided by 4 to fit on the figure.) 
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Figure || displays p(C r (t, t w )) for several choices of the 
ratio t/t w . Interestingly enough, all the curves have a 
noticeable peak at a value of C r that is independent of 
t and t w , with a height that decreases significantly with 
increasing ratio t/t w . The form derived above for tp r (t) — 
<p r (t w ) explains the approximate collapse of p(C r (t, t w )) 
for a fixed ratio t/t w , as shown in Fig. ||for a small value 
of V . (Due to mixing with the stationary part, the t/t w 
scaling worsens when V increases.) Barely noticeable in 
Fig. U is a slow drift of the curves for increasing values 
of t w at fixed ratio t/t w such that the height of the peak 
decreases while the area below the tail at lower values of 
C r increases. This trend leads to the "sub-aging" scaling 
observed for bulk quantities (||| ■ 

Furthermore, the above expression for ip r (t) — f r {tw) 
implies that the PDFs for all of the 28 pairs of times (t, t w ) 
should collapse by rescaling with two parameters: In Ct vp 
and s, corresponding respectively to the nonrandom part 
in ip r (t) — ip r (t w ) and to the width of the random part 
(see Fig. ||). The scaling curve itself gives the PDF for 
X r (t, t w ). The rather good collapse of the curves should 
be improved by further knowledge of C SP . 
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FIG. 4. 
V = 13 3 
n = 1, . . . 



Scaling of p(C r (t,t w )). L = 32, T = 0.72 T c , 
and 28 pairs (t,t w ), with t = 10 4 (2 n_1 ) MCs, 
. 8 and t w = 10 4 (2 11 ™- 1 ) MCs, n w = 1,. . . ,n - 1. 
Ctyp = exp(lnC), where (•••) denotes an average over the 
distribution and s rescales the maximum value of each curve. 
With this choice of V we eliminate the C r < 0, cfr. Fig. ^. 

A local separation of time-scales leads to a 
reparametrization invariant action pd| ] with a soft mode 
that controls the aging dynamics. Our framework, based 
on the analogy with Heiscnberg magnets, predicts the 
existence of a local relationship between C r and Xr, ex- 
pressed by p(C r (t,t w ),Xr(t,t w )) being sharply concen- 
trated along the global x(C) in the (C, x) plane. Under 
additional assumptions, we obtained the scaling behavior 
of p(C r (t,t w )) for all (t,t w ). Our simulations both con- 
firm these predictions and uncover striking regularities 
in the geometry of local fluctuations (l7) . These results 
open the way to a systematic study of local dynamic fluc- 
tuations in glassy systems and they suggest a number of 



exciting avenues for future research. On the theoretical 
side, this framework could be applied to a large variety 
of glassy models, including those without explicit disor- 
der. More interesting still are the experimental tests sug- 
gested by this work. For instance, the local correlations 
of colloidal glasses are accessible experimentally with the 
confocal microscopy technique jq]. Similarly, cantilever 
measurements of noise spectra H allow probing of local 
fluctuations in the glassy phase of polymer melts. These 
are just two examples: any experiment that measures lo- 
cal fluctuations in glassy systems is a potential candidate 
for testing our ideas. 
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